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CN ' We consider Lifshitz field theories with a dynamical critical exponent z equal 

I to the dimension of space d and with a large group of base space symmetries, con- 

cretely space coordinate transformations with unit determinant ("Special Diffeomor- 
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D . phisms" ) . The field configurations of the theories considered may have the topology 

of skyrmions, vortices or monopoles, although we focus our detailed investigations on 
skyrmions. The resulting Lifshitz field theories have a BPS bound and exact soliton 
solutions saturating the bound, as well as time-dependent topological Q-ball solu- 
tions. Finally, we investigate the U(l) gauged versions of the Lifshitz field theories 
p |. coupled to a Chern-Simons gauge field, where the BPS bound and soliton solutions 

pH , saturating the bound continue to exist. 

Csl ' PACS numbers: ll.30.Pb, 11.27.+d 
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CN I of anisotropic scaling of condensed matter physics 2| in other areas of theoretical physics 



X 



rom field theory [3| to gravitation, where it opened a new approach to quantum gravit 
^, and string theory (ADS/CFT) [s^, and back to modern condensed matter theory 
The idea offers an explicit simple realization of systems having different symmetries and 
energy behaviours at short and long distances. Specifically, for the case of quantum field 
theories, an anisotropic scaling of space and time offers the possibility to improve perturba- 
tive renormalizability. These intense activities have achieved many results like a new control 
of Lorentz violations, IR limits of UV renormalizable gravity theories or interesting ABJM 
Chern-Simons duals, but also stumbled on problems, like detailed balance, stability of black 
holes, difficulties with symmetric astrophysical solutions, or the persistence of unwanted de- 
grees of freedom [3]. These can be circumvented for the moment with special readjustments 
and probably are a reflection of general deep problems j^. Therefore, a better understanding 
of the basic field theory input, preferably with exact results even with simplified models, may 
be useful. There have been in fact recent works ^ dealing with solvable (numerically for the 
moment) 3d Lifshitz models with scalars, with gauge couplings or including a Chern-Simon 
term to avoid higher order terms. 
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In this spirit we present here exact solutions and BPS stabihty discussions of specific 
Lifshitz field theories (BPS solitons in Lifshitz theories have also been found in lO]), related 
naturally to special solvable topological soliton models without quadratic terms (quenched) 
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181], with special diffeomorphism (SDiff) invariance, which, as we will show, contain 



in a way some of the anisotropic scaling features, like symmetries. Here the emphasis is 
on nonperturbative solutions carrying a topological charge, which may be of the skyrmion 
or vortex type, and which may be compact or noncompact, depending on the potentials 
considered and their vacuum structure. It is interesting that one obtains exact solutions in 
the specific Lifshitz models considered here. In fact, any exact result in this context can be 
of interest. 

The paper is organized as follows. In Section II we discuss general scaling properties of 
Lifshitz field theories and explain how the topological terms which typically exist in soliton 
models may be naturally incorporated into Lifshitz field theories with dynamical critical 
exponent z equal to the space dimension d. We further show that these Lifshitz soliton 
models contain submodels with invariance under the infinite- dimensional group of area or 
volume preserving diffeomorphisms (SDiff symmetry, in general). In Section III, we study 
in detail a Lifshitz baby Skyrme model which has a BPS bound and both static soliton and 
time-dependent topological Q-ball solutions. In Section IV we study the gauged version of 
this Lifshitz baby Skyrme model with a Chern-Simons interaction of the gauge field. The 
model, again, has a BPS bound and soliton solutions saturating the bound, and this BPS 
bound requires the introduction of a superpotential like in supergravity [l9|, 20 1, analogously 
to the case of the gauged BPS baby Skyrme model Finally, Section V contains our 

conclusions. 



II. LIFSHITZ SCALING AND SDIFF SYMMETRIES 

It is well-known that the presence of higher derivatives in the terms of a field theory 
lagrangian which contribute to the propagator (i.e., terms quadratic in the fields and their 
derivatives) improves the UV behaviour of the theory, because the propagator in momen- 
tum space will contain higher powers of momenta in the denominator. A Lorentz-invariant 
implementation of this idea, however, faces the problem that higher than first time deriva- 
tives in the lagrangian in general spoil unitarity. This lead to the recent proposal jsl, 0] 
of field theories with an anisotropic scaling between space and time coordinates, sacrificing 
thereby (at least in the deep UV region) Lorentz invariance in favor of renormalizability and 
unitarity. Let us briefiy discuss the simple example of a real scalar field in d + 1 space-time 
dimensions. For the moment, we consider the action consisting of the following two terms 
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where is the time derivative of (p, A is the laplacian and A; is a positive integer. For k — 1, 
we recover the standard kinetic term for a scalar field upon a partial integration, but here 
we shall consider the case A; > 1. The action may (and, in general, will) contain more terms, 
but the two terms above are special in that they determine the scaling symmetry and the 
perturbative renormalizability properties of the theory. Concretely, let us discuss the scaling 
behaviour. If we introduce the scaling transformations 

t X K~^x, $ k'^^ (2) 

and require invariance of both terms in the action, then we easily find the following relations, 

2(p — d — k, z — k. (3) 

The two terms $^ and A)'^$ in the lagrangian density corresponding to the above action 
both have scaling dimension d + k. We may add further terms to the lagrangian density, 
but they should not exceed this scaling dimension. Concretely, terms with scaling dimension 
d+k correspond to rcnormalizablc, scale invariant interactions, whereas terms with a smaller 
scaling dimension correspond to supcr-rcnormalizable interactions. 

In a next step, let us make some further assumptions which arc adequate for our purposes. 
First of all, we assume from now on that k = d. This has the consequence that the scalar 
field is now dimensionless, (p — 0. We may, therefore, add arbitrary potential terms V{^) to 
the lagrangian density without spoiling renormalizability. We may also multiply derivative 
terms by functions of $ like, e.g., /($)(— A)'^'$ where, obviously, k' < k must hold. But this 
last observation opens the possibility to include nonlinear terms into the interaction which 
typically show up in field theories with topological solitons like, e.g., the Skyrme model and 
related theories. We simply have to consider more than one scalar field and to act with 
derivatives on them as required by these Skyrme terms. Concretely, we want to consider 
terms which are related to topological charge or winding number densities, which in d space 
dimensions have the typical form 

= /($")e^i-^<^e"i-"'*o'j-,$"i . . . dj^^"''. (4) 

Here, is a function of the fields which depends on the geometry and topology of the 

target space manifold. Further, we assumed that the number of fields is equal to the 
space dimension d, a = l,...,d, which is the most relevant case for topological solitons. 
The square of Q/ is obviously positive semi-definite. In addition, Q'j has scaling dimension 
2d (where we assume k = d) and corresponds, therefore, to a renormalizable, scale invariant 
interaction. 

Up to now, we mainly discussed the effect which Lifshitz-type field theories with 
anisotropic scaling have on perturbative renormalizability, but now we shall slightly change 
our point of view. That is to say, we want to consider Lifshitz field theories not just as 
a fix for renormalization problems, but instead as proper field theories in their own right. 
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which may display new symmetries or new nonperturbative features not present in standard, 
Lorentz-invariant field theories. If the proposal of a field theory with anisotropic scaling in 
the UV is to be taken seriously, then such an investigation is certainly required. Concretely, 
we want to investigate the possible existence of BPS bounds and soliton and Q-ball solutions 
for a class of Lifshitz-type submodels which are characterized by a large group of symmetries. 
The action for the submodels has the generic form 

S=lj d'xdt (^gabm^"^' - X'Q} - /iV($'^)) . (5) 

Here, A and n are real coupling constants, and Qab is a matrix of functions of ("target 
space metric") which obeys det Qab > 0. The two terms gab{^'^)^°'^^ and have scaling 
dimensions 2d, whereas the potential is dimensionless. The absence of a term of the 

type $(— A)*^^ implies that these submodels are certainly not useful for the discussion of 
perturbative features, but they may be useful to understand some nonperturbative prop- 
erties of Lifshitz-type theories. The reason to choose the above actions is that they have 
an infinite-dimensional group of symmetries. They are invariant under base space SDiff 
("Special Diffeomorphism" ) transformations, i.e., diffeomorphisms of the space coordinates 

— )■ x'^ = y^{x^) which leave the volume form d'^x invariant (i^e., diffeomorphisms with 
unit determinant). We remark that in Horava-Lifshitz gravity [4] the action is invariant 
under foliation-preserving diffeomorphisms, therefore our above actions may be seen as toy 
models for the study of certain nonperturbative issues of that proposal after a partial gauge 
fixin g (s etting the lapse function equal to one and the shift functions equal to zero, see 
also 21j, 221 for further discussion). The infinitely many symmetries (and corresponding 
conservation laws) of the above actions ([5]) imply that these theories are integrable in the 
sense of generalized integrability 123], and, therefore, one expects the existence of exact 
nonperturbative solutions. 

If we restrict to static field configurations, then the corresponding energy functional 

^[0-] = 1 y d^xdt (A'QJ + /iV($'^)) (6) 

has even more symmetries. Indeed, the first term is invariant under the group of target space 
SDiff symmetries, that is, under the target space diffeomorphisms which leave the target 
space volume form = /($")(i$^ . . . d<^'^ invariant. The potential, in general, breaks some 
of these SDiff symmetries, but for sufficiently symmetric potentials an infinite-dimensional 
symmetry group remains. For potentials, e.g., which only depend on the modulus all 
SDiff transformation which do not change this modulus remain symmetries of the full static 
energy. It is interesting to compare the symmetries of the actions ([5]) with the symmetries 
of the BPS-Skyrme type theories - 18] whose actions are simply the Lorentz-invariant 
generalizations to time-dependent field configurations of the above energy ([6]), i.e., their 
actions have the form 

^ = 1 y df^xdi {-\^{Q^ff - /iV(<l>")) (7) 
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where 

is a Lorentz-covariant topological current. The static energy functionals of the two types 
of theories are the same, therefore also their symmetries coincide. For the time dependent 
situation, on the other hand, the roles of the two groups of symmetries are inverted for the 
two types of theories. That is to say, while the actions still have the base space SDiffs 
as symmetries (these are, therefore, Noether symmetries), but no longer the target space 
SDiffs (more precisely, the subgroup which leaves the potential invariant), the situation is 
exactly the other way round for the theories ([7]). These actions are still invariant under the 
target space SDiff subgroup, but no longer under the base space SDiffs. Another difference 
between the Lifshitz-type theories (I5j) and the BPS-Skyrme type theories is that the 
former have a well-defined Cauchy problem, whereas the latter ones don't (i.e., for some 
initial data their time evolution is not well-defined, see e.g. [l3|, |18|). 

Finally, we shall restrict to the case of two space dimensions, d = 2, for the concrete 
examples considered below. The main reason is that calculations are simpler there, although 
for most results the generalization to higher dimensions should pose no problem. Besides, the 
field configurations in the examples considered will have the topology of (baby) skyrmions. 
We could also easily generalize to fields with the topology of vortices (or monopoles in 
d = 3 dimensions) - indeed, for the BPS Skyrme type models there even exists an 
exact topological duality between theories and solutions with Skyrme (soliton) topology and 
theories and solutions with vortex (or monopole) topology [l8^. But for reasons of simplicity 
we restrict our examples to baby skyrmions. 



III. THE LIFSHITZ BPS BABY SKYRME MODEL 



For baby skyrmions 2J], the target space manifold is the unit two-sphere S"^. This implies 
that the fields a = 1, 2 should be interpreted as the real and imaginary part, u = <l>^+i$^, 
of a stereographic projection (or CP(1)) coordinate u E C. The target space metric and 
area two-form are, therefore, 

4 4 
and the topological current is 

and provides the degree of the map defined by the Skyrme field $"(a;) (topological charge 
or winding number) according to 

deg[0] = / d^xq e Z (11) 
47r ./ 
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where we use the short-hand notation q = q^. In the following, however, we prefer to 
parametrize the target space two-sphere by a three-component unit vector 0, i.e., 0^ = 1 
related to u = + by stereographic projection. In this parametrization, there are no 
nontrivial target space metric and area factors, and the nontrivial target space geometry is 
taken into account by the constraint 0^ = 1. E.g., the topological current reads 

e^i^A?^ = = 0- (<9^0 X (9^0) . (12) 



A. The model 

In contrast to the BPS baby Skyrme model [3] - [isl, its Lifshitz version has the standard 
time derivative part in the Lagrangian (but remember the constraint 0^ = 1), while the term 
containing the space derivatives remains unchanged, 

C = j{do^f-^q'-fi'V{n-$). (13) 

Here, the potential V{n ■ 0) breaks the S0(3) target space symmetry (rotations of 0) down 
to an S0(2) symmetry (rotations about the n axis). Further, the potential obeys V > 0, 
V{1) = 0, and n is the vacuum vector which must approach in the large distance limit for 
finite energy field configurations, lim|^|_j.oo = n. Further, 7, A and /i are real constants. 

In the static limit, the model is identical to the BPS baby Skyrme model and, there- 
fore, possesses the same exact static solutions, which, depending on the potential V, are 
compactons or usual infinitely extended solitons with power-like or exponential localization. 
Let us briefly summarize some known results [l5|. For the family of old baby potentials 
(one- vacuum potentials) 

V-C-^X (14) 



2 

where a G (0,oo), we get compactons (a G (0,2)), one exponential solution (a = 2) and 
power-like localized solitons (a > 2). We remark that the effect of equivalent potentials for 



generalizations with a CP^ target space have been studied recently in 25j. On the other 
hand, for the family of the new baby potentials (two-vacuum potentials) 

v={i-{<j>rr (15) 



one flnds only compactons for a G (0,2). For other potentials [26|, [27| and exact solutions 



see 
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While the static spectra of the BPS baby Skyrme model and the Lifshitz version coincide, 
their time dependent solutions as well as the dynamical properties of the static configurations 
differ dramatically. In the case of the BPS baby Skyrme model, the time evolution is a highly 
non-trivial problem since the equations are not globally hyperbolic. Hence, for some initial 
configurations time evolution may be not globally well defined and could lead to singularities. 
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However, due to the fact that the Lagrangian depends only on first time derivatives squared 
one can still define a proper Hamiltonian and (which is important from the point of view of 
the application to the nuclear physics) perform the semiclassical quantization. Obviously, 
the Lifshitz version has a well defined dynamics which can be at least treated using numerical 
methods. 



B. Dynamics 

The full time dependent equations of motion read 

- 270 X dl$ + X $V' + X^ViKi = (16) 

where 

i?, = V,0(0-(V,0x V,0)). (17) 

Let us observe that the Lagrangian looks like the classical membrane Lagrangian with a 
nontrivial potential 0] 

£ = ^{doX'^f - X^}2 - /iV(X'^) (18) 

where 

{X^X''} = e,,.9,X"9,X'' (19) 

is the Poisson bracket on the base space while X"- = (f)"-. Then, the static Bogomolny 
equation is just 

^abc^j^a^ ^ TVX" = (20) 

and is identical to the BPS equation of the BPS baby Skyrme model. On the other hand, 
for the usual membrane (with potential being in fact a constraint - see, e.g., V = (X")^), 
the second order dynamics can be reduced to a first order Nahm equation 

doX" = e"''"{X^ X^} + X". (21) 

Due to the similarity with our Bogomolny equation, one might hope that also here the 
dynamics could be integrated to a first order "nonlinear" Nahm equation 

9oX" = e"^"{X^ X^} + WX^. (22) 

This is, however, not true, and the djTiamics of the two theories is completely different. The 
easiest way to see it is to multiply the last formula by X", i.e., 

X^doX" = e''^''X''{X\ X"} + W. (23) 

If we now assume that the X"^ span a two-sphere, then X^X"^ = 0, the r.h.s. is exactly the 
Bogomolny equation and, as X'^(e"^'=X"{X^ X^}) = e'^'"={X^X^} we find that 

e^'"{X^, X''} + WX'' = (24) 
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and, hence 

doX" = (25) 

i.e., no dynamics in the system. This just corresponds to the well-known fact that the sphere 
is a vacuum solution of the membrane action. On the other hand, we expect nontrivial 
dynamics for our Lifshitz model with the target space manifold S^, so the dynamics (time 
dependence) of the two theories is certainly different. 

C. Spinning compacton solutions 

It is possible to find exact time dependent solutions assuming the following axially sym- 
metric ansatz {x^ — r cos (f), x"^ — r sin (p) 



<P{r,(l),t) 



^ sin /(r) cos(n</) — cut) ^ 

sin /(r) sin(n0 — ut) 
\ cos/(r) J 



(26) 



where u is a real parameter of the internal rotation. Here, the integer n is equal to the 
topological charge of the skyrmion provided that the profile function /(r) satisfied the cor- 
responding skyrmion boundary conditions. The resulting Lagrangian reads 

L^2t: J dx {4jLJ^h{l -h)- 2\^n^hl - /iV(/i)) (27) 
where we have assumed fi = (0, 0, 1) =^ n • $ — (f)z and used the new target space variable 

h^\{l-h) (28) 
Moreover, x — /2. The pertinent equation of motion is 

-4AV/i^^ + /xV^-47a;2(l-2/i) = (29) 
where the nontrivial topology of the Skyrme field enforces the following boundary conditions 

h{x = 0) = 1, h{x = X) = 0, h^{x = X) = (30) 

where X can be finite (compactons) or infinite (usual soliton). Taking into account the 
boundary conditions, this equation can be easily integrated to 

2AV/i2 =yuV- 476^2/1(1-/1). (31) 

Of course, one can integrate this formula again and easily express the solution in terms of 
the integral depending on the potential V . Hence, similarly to the static case we find exact 
solutions. 



9 



It should be stressed that in order to have a topologically nontrivial solution h must cover 
the segment [0, 1] which results in a condition for the effective potential. Namely, 

Veff ^^V- 4^h{l - h) (32) 

must be positive definite on the whole segment h e [0,1]. This can restrict the possible 
values of a; as well as exclude the existence of the spinning solutions completely. 

Proposition 1: Generic spinning solutions are of a compact type while the critical 
solutions are compactons or infinitely extended solitons (exponentially or power-like 
localized) . 

Proof: The approach to the vacuum, and therefore the qualitative type of a solitonic 
solution, is controlled by the asymptotical behavior of the effective potential near the 
vacuum. Here, in order to have finite energy configurations, we assume that V (and in 
a consequence Vg//) must possesses one or two vacua located at /i = (respectively at 
= 0, 1). As the second term in the effective potential is fixed, it is enough to analyze the 
Taylor expansion of the original potential at /i = 

y = co/i^ + o(/i^). (33) 

If c < 1 then also the leading term of \4// goes sublinearly with the field. Due to that one 
obtains a compacton (which ceases to exist for sufficiently large uj). Moreover, there are no 
spinning solutions for any potentials with c > 1. In this case the effective potential has the 
linear leading linear term, however, with a negative prefactor. In other words, at /i = the 
effective potential always gets a negative value. The last possibility when c = 1 is slightly 
more subtle. Generically one finds again compactons as is linear in h. However, there 
is a unique value of the frequency ui (we called it the critical frequency ujcrit) for which the 
linear leading term in V cancels with the linear term of the second part of 

ujcrit = (34) 
Then, one has to take into account the next term in the expansion 

V ^coh + ciN^^ +o{h'^^). (35) 



If (ii < 2 then we get a compacton also in the critical limit (ci must be positive). For di > 2 
we find an exponentially localized solution as [A'-yu'^ / n^)h'^ becomes the leading term in Vefj. 
In the case di = 2 there arc two possibilities: if there arc no other terms in the expansion 
we get a soliton with the exponential tail (under the condition that Ci + Ajcu'^^^^/ n'^ > 0). If 
further higher terms in the expansion exist, power-like localized solutions are observed for 
ci = — cq. For CO > LOcrit no solutions exist. Obviously, such non-compact solutions are quite 
special (not generic) since they emerge for a unique value of the spinning frequency and/or 
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for very special potentials. 

Observation 1: The old baby type potentials generating static BPS baby skyrmions with 
exponential or power-like localization do not lead to solitonic spinning solutions 

As an example with exact solution let us consider the old baby potential V = h. 
Then we find 



hi = ^h{a + f3h) (36) 

77/ 



where 



«^ = ^(/^^-47-^), = ^ (37) 

In order to have a topologically nontrivial solution h must cover the segment [0, 1] which 
requires a > 0. Hence we get an upper bound for the frequency of the oscillation 

< i^crit = ^ (38) 

Then the exact compacton solution is given by 

(39) 

where the compacton boundary 

p a 

In the critical point where uj = cOcrit and a = we find a non-compacton exponentially 
localized solution 

hcrit = e-^"" (40) 

This is consistent with the fact that the size of the compacton grows to infinity as we 
approach the critical frequency. It is an interesting observation that the localization type of 
the solution depends not only on the original potential but also on the value of the oscillation 
frequency as both contribute to the effective potential. 
Similarly, for the new baby potential V = h{l — h) we find 



hl-J^. - 47^') Ml - h) (41) 

which represents a compacton with the same bound on the frequency. The corresponding 
compact baby Skyrmion reads 

f sin2^(x-X) x<X 

I x>X ^ ' 

and 

a 
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Using the first order profile equation we can rewrite the total energy integral of the spinning 
solutions as 

E = 2Ti j dx (47u;2/i(l - h) + 2\^n^hl + /iV(/i)) = 4^ j dxV{h) (43) 
which for the above considered old baby potential gives 




Bold = —7= I A / — A / IH arc sinh 

while the new baby potential leads to 

Er, 




-'new 



2a a/2 li 47a; 



Hence, both case lead to linear energy-topological charge relations, with the critical fre- 
quency fixed by the parameters of the model (and independent of the topological charge 
Q = n). This should be contrasted with the baby Skyrme model and its energies of spin- 



ning solutions [27[, where the energy-charge relation is linear only approximately. In fact, 
E/Q depends on the charge which leads to some instabilities for axially symmetric higher 
charge solutions. Moreover, there has been observed a symmetry breaking transition from 
the axially symmetrical baby skyrmions to discrete symmetry ones if the values of the model 
parameters are varied. The linear energy-charge relation found here may point toward the 
stability of the axially symmetric solutions. 

As we see, the Lifshitz BPS baby Skyrme model is in a sense a non-trivial combination of the 
usual Skyrme model and its BPS version sharing some properties with former while others 
with the latter model. 

In common with the BPS baby Skyrme model, we find that the spinning solutions are axi- 
ally symmetric with energy proportional to the topological charge. There is also a maximal 
critical frequency above which solutions cease to exist. However, in the BPS baby Skyrme 
case we still have compacton solutions while here the size grows to infinity. Another dif- 
ference is that here Ucrit is topological charge independent, while in the BPS baby model 
it goes as ^/n. Further, in the BPS model the size of the compactons in fact becomes 
smaller and smaller as we increase the frequency. Hence, in this aspect the Lifshitz model is 
more like the usual baby Skyrme model. As a general conclusion, the different structure of 
(time-dependent) Lifshitz and Lorentz-symmetric field theories may lead to a qualitatively 
different behaviour of their time dependent solutions (for recent results on time- dependent 



solutions of Lorentz- violating nonlinear field theories we refer, e.g., to 29| ) 



IV. THE GAUGED LIFSHITZ BPS BABY SKYRME MODEL 

The potential V{n ■ 0) leaves an S0(2) ~ U(l) subgroup of the target space S0(3) 
symmetry intact, therefore a natural way to couple the Skjnrme field to an electromagnetic 
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gauge field is by gauging this residual U(l) symmetry. We shall require that the introduction 
of the gauge coupling via minimal substitution does not spoil the scaling properties of the 
theory, which implies that the gauge fields must scale under the anisotropic scaling like the 
corresponding partial derivatives. That is to say, Aq should have scaling dimension z (where, 
in our case, z = d = 2), whereas the spatial components Aj must have scaling dimensions 
equal to one. But this results in a Maxwell term where different contributions have 
different scaling dimensions. Specifically, terms with time derivatives do not correspond to 
renormalizable interactions. In d = 2 space dimensions, one way to avoid this problem is 
by introducing a Chern-Simons term instead of the Maxwell term. It may be checked easily 
that all contributions to the Chern-Simons term have scahng dimension z + d — 2d — 4: and, 
therefore, correspond to a scale- invariant, renormalizable interaction. 

A. The model coupled to Chern-Simons electromagnetism 

The model we want to consider is the Lifshitz BPS baby Skyrme model of the last section, 
minimally coupled to Chern-Simons electromagentism. The corresponding lagrangian is 

C = 7(^00)' - y + ^e'^'^'^F.^A^ - i,'V{n ■ 0) (44) 

where 

Q^^ = • {D^(j) X D^(p) = q^^ + A^d,,{n ■ $) - A^d^{n ■ $) = e^^A*?^, (45) 

i.e., 

Qi. ^ le^-Ag^^ ^ + e^^'^A^d^ifi ■ $), (46) 

and Q = Q'^. Here, p is a new coupling constant measuring the coupling of the Chern-Simons 
term, and D/j, is the usual covariant derivative 

Dj^di,$+Ai,nx$. (47) 

B. The BPS bound 

Using the short-hand notation Q = Q^, q = q^, the lagrangian density may be re- 
expressed like 

\ 2 

C = ^(^$^ + 2Ao$- (n X 0) + Al (l - (n • 0)^)) - —Q^ + p{AoB + e^jEiAj) - /^V (48) 
where we made the Aq dependence explicit and introduced the electric and magnetic fields 



Ei = doAi - diAo , B = Fi2 = diA2 - 82 Ai. 



(49) 
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We want to study the energy functional for static fields, but before this we have to solve the 
constraint equation (the Gauss law) of the model. The Gauss law is 

^ -5,^ = 27 (0- {n X 0) + Ao (1 - {n-$f))+2pB ^ (50) 

or, for static configurations 

B = - 

P 



i? = -^Ao(l-(n-0)). (51) 



The static lagrangian is 

C = ^Al (1 - (n ■ 0)') - + p{A^B + t,,E,A,) - /iV (52) 

and, using 

eijEiAj = -eij{diAQ)Aj = -O^^ei^AQA^) + AqB, (53) 

skipping the total derivative, and eliminating Aq via the Gauss law, we get the static energy 
density 

£ = -C = P- + ^Q^ + p?V{n ■ $) (54) 

7 1 - (n ■ 0)2 2, 

which is positive semi-definite. The total derivative —di{eijAQAj) which we omitted does 
not influence the Euler-Lagrange equations, but it might produce boundary contributions 
to the static energy, in which case the two expressions fl52l) and fl5^ for the static energies 
/ dPx£ = — J dPxC would give different results. We shall find, however, that BPS solutions 
always obey lim|^|^oo ^0 = 0, so the boundary term does not contribute and the two energy 
expressions provide the same energy. 

Now we assume n = (0, 0, 1) ^ n ■ = 03 and use the new variable h defined in (128|) . 
then the static energy functional is 

E = [d'x( + ^TTi^ + p'V{h)\ . (55) 



2^ 47/1(1-/1) 

This is similar to the energy of the gauged BPS baby Skyrme model 17|, and also the 
BPS bound is similar. Indeed, we introduce a superpotential W{h) and start with the 
non-negative expression 

[^(Q= + wi) + ^ + l^Mi - m-) (56) 

+A' {QWh + 2BW)\ . (57) 

Using Q = q — 2eijAidjh the last line may be written as a topological term plus a total 
derivative, 

QWh + 2BW = WhQ - 2dj{eijAiW). (58) 
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We will find that Ai does not vanish in the limit x — )■ oo, therefore the total derivative can 
be omitted only provided that W obeys lim^_>oo W = 0. But W depends on h only, and h 
must take the vacuum value h = at infinity, therefore the condition which W has to obey 
is 

W{h = 0) = 0. (59) 

We shall find that for the class of potentials V{h) we consider in this paper, the corresponding 
superpotentials indeed always obey the above boundary condition. 

Assuming this, and comparing the non-negative expression with the energy, we find that 
the energy obeys the BPS inequality 

E>\^\ J d^xqWh\ (60) 



provided that the superpotential W obeys the modified superpotential equation 

+ ^h{l - h)W' = /xV (61) 



together with the boundary condition ( 159|) . Finally, the BPS bound is topological because 
q is the pullback of the area two-form dQ on the target space two-sphere S^, and, therefore, 

SxqWh = n dVLWh = n dip d(l)^Wh = 47m / dhWh = AnnWil) (62) 
J Jo J-i Jo 

where n is the winding number (degree) of the map defined by the Skyrme field configuration 
0. 

The BPS bound is saturated for fields which obey the corresponding BPS equations 

Q = -Wh, B = -^K^ - h)W. (63) 

It may be proved that the BPS equations imply the static Euler-Lagrange equations. For 
general static configurations, the proof is rather lengthy and is similar to the proof in 17 . 
Here we shall give the explicit proof for axially symmetric configurations in the next section. 

We remark that, in principle, there exists a further term with Lifshitz scaling dimension 
equal to z + d = 2d = A which we might include, namely the gauge potential coupled (non- 
minimally) to the topological current. The corresponding contribution to the action will 
be gauge invariant only provided that the current is both conserved and gauge invariant, 
so neither q^ (which is not gauge invariant) nor (which is not conserved) may be used. 
There exists, however, a gauge invariant and conserved topological current, namely = 



Q'^ + {l/2)e^'^''F^p{l — n ■ (p) [30|, therefore the term p J d'^xdtQ^A^ is gauge invariant. The 
resulting theory with this term included still has a BPS bound and BPS equations in terms of 
a superpotential, and the structure of the superpotential equation is similar (first order and 



quadratic). Concretely, the superpotential equation has the form AW'"^ + 2BWW' + CW'^ 
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V, where A, B, C are (quite complicated) rational functions of h and define a non-negative 
quadratic form, i.e., AC — B'^ > 0. BPS soliton solutions could then be found, e.g., by 
numerically integrating the resulting three first order equations (the two BPS equations and 
the superpotential equation). For reasons of simplicity, however, we will not further consider 
this general case here and set p = in the rest of the paper. 

C. The axially symmetric ansatz 

Next we assume n = (0, 0, 1) and the standard static ansatz 



^ sin /(r) cosn0 ^ 

sin /(r) sin 720 
V cos/(r) J 



Aq = nb{r), Ar = 0, A^ = na{r). (64) 



Then the magnetic field is 



B = \e'^^F,, = e'^^d.A, = (65) 



The resulting Lagrangian is 



L = 2nn^ f rdr ( 76' sin^ / - ^2 sin^ /4(1 + a)^ + ^(a,6 - Ka) - ^Vif)] . (66) 
J \ ^ ^2 ^ J 

We introduce the variable x = and the new target space variable h defined in fl28l) . 
where ^ ^ 

h = -il-cosf) ^ /i^ = -sin//^, (67) 

L = 2Tm^ J dx ^476^/1(1 -h)- 2X'^hl{l + a)2 + p{a^b - h^a) - ^V{h)^ . (68) 



then 



+ a)2) - -^Vf, + ^b\l -2h) = (69) 



The static field equations read 

4A2n2'" ' X2 

b, = -—hl{l + a), (70) 

and the Gauss law is 

a^ = -^bh{l-h). (71) 
Using the Gauss law to eliminate b we get the positive definite Hamiltonian 



H = -L = 2nn' / dx ( 2X'hl{l + af + ^^^^^ + ^,V{h) ) . (72) 



2 
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In a next step, let us prove that, for the axially symmetric ansatz, the BPS equations 
together with the superpotential equation and the Gauss law imply the static second order 
equations. The BPS equations for the axially symmetric ansatz read 

2nh,il + a) = -Wh (73) 
na^ = ^^(1 - h)W. (74) 

We remark that the superpotential equation fl6ip is invariant under the sign change W — )• 
—W. We chose the sign such that W and Wh will be positive in general. This follows from 
the fact that hx is negative, in general, because h{0) = 1, h{oo) = 0. 

Using the second BPS equation ( 17il) and the axially symmetric Gauss law (17T1) to eliminate 
ttx we easily find 

\ 2 \ 2 \ 2 

b = —W bx = —Whhx = -2-hl{l + a), (75) 

np np p 

i.e., Eq. ( 170]) . where we used the first BPS equation fITH]) in the last step. Finally, the field 
equation (|69|) may be transformed into the h derivative of the superpotential equation by 
the following series of steps. 



2 



3. (ft,(l + af) + ^6^(1 - 2ft) = 

1 \ 2 2 

-— {WhhKil + a) + Wnax) + \^W\l - 2h) = j^V^ 

2n n'^p^ 4n^A^ 



2 p^ 



dh ^Wt + ^h(\ - h)W' = p'Vh (76) 



where the last line is just the h derivative of the superpotential equation, and we used the 
BPS equations and the Gauss law in several instances. Reading the above derivation from 
below to above, it follows that also the first field equation (1^ is a consequence of the BPS 
equations and the superpotential equation. 



D. The condition = 0) = 



Now, let us demonstrate that for the class of potentials we want to consider the corre- 
sponding superpotential which gives rise to the BPS solutions must always obey the bound- 
ary condition W{h = 0) = 0. The class of potentials we permit have their vacuum at /i = 0, 
i.e., V{h = 0) = 0. We do not exclude the possibility of more vacua, but it is always the 
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vacuum ai h = which is approached in the hmit of large \x\, i.e. \im\^\_^oo h{x) = 0. 
Further, we assume that the vacuum at /i = is approached hke 

V{h) ~ /i" , a > 1. (77) 

For a > 1, the condition W{0) = is a direct consequence of the superpotential equation 
( I6T1) . Indeed, the superpotential equation near h = implies that Wh ~ h^°'^'^\ which a 
priori seems to allow for the superpotential W ~ -)_ where Wq is a constant. But 

let us study the h derivative of the superpotential equation near h = 0, 

X^WhWhh + ^((1 - ^h)W^ + 2h{l - h)WWh) = /i'a/i"-^ (78) 

For a > 1, both the r.h.s. of the above equation and the product WhWhh are zero at = 0, 
which implies W{h = 0) = 0, which is what we wanted to prove. 

The case a = 1, on the other hand, is more involved. It is, at the same time, important, 
because it includes well-known potentials like the "old" and "new" baby Skyrme potentials. 
In the above equation (|78l) . both WhWhh and h°'~^ give finite, nonzero contributions at 
h = for a = 1, and the superpotential equation allows, in fact, a one-parameter family of 
solutions parametrized by the boundary value Wq = W{0). We shall demonstrate, however, 
that only the case Wq = leads to a genuine soliton solution with finite energy, whereas 
other choices lead to formal solutions with singularities or infinite energy. For simplicity, we 
will restrict to the case of spherically symmetric field configurations h{x), a{x), i.e., to the 
variational problem defined by the energy functional (|72l) . It turns out that the discussion 
is simpler directly for the second order system resulting from this energy, where the relation 
to the first order (BPS) system will be obvious. The field equations resulting from (172|) are 
(equivalently they may be derived from Eqs. (l69l) . (170!) and (17T]) by eliminating b with the 
help of the Gauss law) 

4A^4(Ml + ««+y^f^|^«I = ^V.. (79) 

,XVi(l^a)^^8,[^). (80) 

As a next step, we need the information that for V ~ /i the solitons are compactons which 
approach their vacuum values at a finite distance xq. Concretely, the approach is quadratic, 
i.e., like h ~ h2{xo — x)^. As a consequence, the boundary conditions for the variational 
problem are 

h{0) = 1, a(0) = 0, h{xo) = 0, h'{xo) = 0, a{xo) = 0. (81) 

On the other hand, we have a total number of five free constants (four integration constants 
and the compacton radius Xq), so free constants and boundary conditions seem to match 
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exactly, leading to at most one solution. This conclusion is, however, not correct. The 
problem is that the condition a'(xo) = is not an independent condition but, instead, a 
consequence of the field equations and the remaining boundary conditions. So apparently 
we found a one-parameter family of solutions (five free parameters and four independent 
boundary conditions). The problem with this family of apparent solutions is that different 
solutions will, in general, have different energies. Indeed, for a fixed value of the compacton 
radius xq the number of free parameters and boundary conditions match, so we expect 
at most one solution. Different values of Xq, however, correspond to different variational 
problems (the same equations and boundary conditions, but on different intervals), so there 
is no reason why the energies should be the same. This would imply that the original 
variational problem is ill-defined. The resolution of the puzzle is that we have to impose 
one more condition for a genuine compacton. A genuine compacton is not a solution on the 
finite interval [0,Xo]. Instead, it is a solution on the full half-line which allows to connect a 
nontrivial solution for x G [0, xq] with the vacuum solution for x > Xo- The most direct way 
to derive this additional condition is by inserting the power series expansion 

oo 

h = J2hk{xo-x)'' (82) 

k=2 

oo 

a = J2(^kixo-x)'' (83) 

fc=0 

into the field equations fl79|) . flHOl) . As we are expanding near the vacuum h = 0, we replace 
the potential V hj V ^ h and \4 by \4 — )■ 1. We find that ao remains undetermined, 
whereas ai = and 02 = 0. The next condition is 



%l + hi + 8X\1 + aofh^] = 0. (84) 

47 \ / 



This condition should be interpreted as a determining equation for /i2 for given values of ao 
and 03. It is a cubic equation in h2, so it will, in general, lead to three solutions (roots) for 
h2- It will, therefore, allow to join different solutions at Xq, corresponding to the choice of 
different roots h2. The vacuum solution /12 = 0, however, may be found among the roots only 
for 03 = 0. For a compacton solution we must, therefore, require the additional condition 
03 = 0. For this choice we may join the nontrivial solution 

(85) 



8n2A2(l + ao)2 



for X < Xq with the vacuum solution /12 = for x > xq- In addition, we have now two 
free constants (oq and xq) to satisfy the two remaining boundary conditions h{0) = 1 and 
a(0) = at the center x = 0, so we expect at most one solution, and, therefore, a well-defined 
variational problem. 
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Finally, it easily follows from the BPS equation Eq. (17^ that 03 = implies W{h = 0) = 
which, as a consequence, continues to hold also for potentials which are linear in h near 
the vacuum h = 0. It further follows from Eq. ( 175|) that W{h = 0) = implies b{xo) = 
(or lim^.__>oo K^) = non-compact solitons) and, therefore, lim|;j;|^oo ^oi^) =0. As a 
result, the total derivative terms in Eqs. f l53|) and f l58|l do not contribute to the energy, as 
announced. 



E. The magnetic flux 

It is possible to find an exact expression (i.e., in terms of the superpotential W) for the 
magnetic flux 

$ = y rdrd(j)B = 27m J dxttx = 27ma{xo) = 27rnaoo (86) 

where Xq is the position of the boundary of the soliton which can be finite in the case of 
compactons as well as infinite for the usual exponentially or power-like localized solutions. 
From the BPS equations for the axially symmetric ansatz we get 

_ 87A^/i(l - h )W.^ 

which can be integrated to 



-hx, (87) 



dx ln(l + a) = ^-^dxF{h) (88) 
where 



It gives 



lnC(l + a) = ^F(/i(x)) 1 + a(x) = C'-^e^^^''^"" (90) 

P 



where the integration constant C can be eliminated using the boundary conditions h{x 
0) = 1 and a{x = 0) = 0. Then, 



Hence finally 



lnC=^F(l) =, C = e'f'^'^ (91) 



-1 + e^^^'^ (92) 



where h{xo) = and F{h = 0) = have been used. As we see, the final expression is quite 
similar to the case of the gauged BPS baby case with the Maxwell gauge field. However, 
the specific form of the function F{h) is different which could potentially lead to different 
properties as, e.g., quantization of the magnetic flux. 

The quantization of the flux would appear if the function F tends to infinity at h = 1 which 
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obviously requires that the integrand has a singularity (1 — h)~^ or stronger. The simplest 
realization of this case is given by the superpotential W having the following properties: 

(i) W{0) = and no other zeros for W 

(ii) Wh^{l- /^)^ where a > 2 as /i -> 1. 

As an example one may consider the superpotential W in the following form 

W = + _ = /iM + _ 93) 

3 2 5 V3 2 5 y ^ ^ 

which has only one zero for /i = on the segment h G [0, 1]. Then 

Wh = h\l - hf. (94) 

From the superpotential equation we may easily find the potential V corresponding to this 
choice. It is a well behaving potential with two vacua located at the ends of the segment. 
Now 

which is singular at h = 1. It means that -F(l) = oo and 

aoo = -1. (96) 

So, at least for this potential V we get quantization of the magnetic flux for all values of the 
coupling constants. However, one must check whether this superpotential leads to solitonic 
solutions with the prescribed boundary conditions. We will see that this is not the case. 
It can be easily shown from the first BPS equation (173|) . In the vicinity of the second vacuum 
/i = 1 (at X = 0) the gauge field vanishes and we arrive at the simple equation 

2nh^ ~ (1 - hy-^ (97) 

where we used the superpotential given above. The pertinent solution at the origin reads 

1 — h —— — -r. 

Hence, h = 1 can be reached only for x oo, which is incompatible with the boundary 
condition at the origin. In other words, the only solution with our boundary conditions is 
the trivial (vacuum) solution. Therefore, the fiux cannot be quantized. We remark that in 
this derivation we only needed the behaviour of Wh near h = 1, which, in turn, was necessary 
for the condition -F(l) = oo ^ a^o = —1, so the impossibility of fiux quantization seems to 
be a general result. 

As usual for the Chern-Simon type models, the magnetic fiux is related to the electric charge. 
To see it, let us write the field equation for the gauge field as 



(98) 
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where the current j'^ contains only the gauge field and not its derivatives. Obviously the 
current is conserved = 0. Knowing that 



^--t/' = '''' 



we get 



q> = I (i^xB = -j- J d^xf = (100) 



where Q is the electric charge. 



F. Gradient flow 

Here we want to demonstrate that the BPS equations for the axially symmetric ansatz 
may be expressed as a gradient flow equation w.r.t. a certain target space metric, exactly like 
in the cases where the superpotential and the corresponding BPS equations are derived via 



"fake supergravity" for self-gravitating domain wa' 



Is, inflaton scalar flelds, or extremal black 



20|. Let us start with the BPS equations 



holes with some (hyper-)spherical symmetry 19| 
in the axially symmetric ansatz 

2nK{l + a) = -Wh (101) 
47A2 

na^ = ' h(l-h)W (102) 

P 

Now, we introduce a new target space variable a + 1 = and a new superpotential 
W = e^W{h). Then we arrive at 

hx = (103) 

= - h)Wu (104) 

One can write such BPS equations as flrst order gradient flow equations 

r - G^'^ = (105) 

where = {U,h) are the target space coordinates. The "effective metric" on the target 
space is 







Gat = '''' f -'^^ , . (106) 
\ —2n J 

As in the case of the BPS baby Skyrme model coupled to the Maxwell gauge fleld, we 
flnd a Minskowskian signature of the gradient flow target space. This resembles the case 
of extremal black holes but is different from the cases of domain walls and inflaton flelds, 
which lead to an Euclidean signature. However, our metric is not constant. It becomes 
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singular at the north and south poles of the S*^ 3 (j). We remark that a detailed discussion 
of g radient flows and the related supertpotentials in Lifshitz field theories can be found in 
22l | where, however, the flow is parametrized by euclidean time and not by a space variable 
(like X = in our case). 



G. Exactly solvable superpotential equation 



For the new baby potential V = h{l — h) we can solve the superpotential equation 
completely. We have 

—Wh + - h)W^ = fi^h{l - h) (107) 



^Wl = h{l - h)^? {l - i^iy^) (108) 

A dW , 

±^h{\-h)dh (109) 



V2/^yrr|s^ 



Hence, 



arcsin — W = ±- 



2(1 - 2h)^/h{l -h) + arcsin(l - 2h) + C (110) 



or 

W 



h l^f sinf [2(1 - 2h)Jh{l -h) + arcsin(l - 2h) + c]] (111) 

2A^y7 \2V2p L jy 

where C is an integration constant fixing different W{0) values. We already know that 
W{0) = is the only acceptable value, which implies C = —{n/2). 

From the explicit expression above it is obvious that a superpotential with the right 
boundary condition W{0) = exists for arbitrary values of the coupling constants. This 
does not imply, however, that BPS soliton solutions will exist for all values of the couplings. 
Indeed, for sufficiently large values of ^5^, the derivative Wh will have a zero in the interval 
h G (0, 1), and in this case a BPS soliton does not exist. The argument goes like follows. Let 
us assume that Wh has exactly one zero in the open interval h G (0, 1). Then Wh is either 
positive below the zero and negative above, or vice versa. In any case, it has different signs 
near h = and near h = 1. But via the BPS equation 2nhx{l + a) = —Wh this implies that 

has different signs near a; = and near x = 00, which is incompatible with the conditions 
h{0) = 1, h{oo) = and h G [0, 1]. In this argument we used that 1 + a{x) > V x, as 
follows from Eq. fl90|) . This argument does not exclude the possiblity of BPS solitons for 
superpotentials such that Wh has an even number of zeros in the interval (0, 1). But for the 
exact superpotential of this section, this last possibility can be excluded, too. In order to 
see it, we may consider the expression for the total flux which is completely fixed by the 
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value of the previously introduced function F{h) at h = 1 (see Section IV.D). In the case of 
the exact superpotential we get 
1 

F(l) = — / dh^h{l - h) tan (k [2(1 - 2h)y/h{l - h) + arcsin(l - 2h) - -1) (112) 
8k J Q V L 2J/ 

1 

= — dh^Jh{l - h) tan{Kp{h)) (113) 
8fi Jo 

where k = X^/{2y/2p) and p{h) = 2(1 - 2h)^Jh{l - h) + arcsin(l - 2h) - f . Surprisingly 
the integral is analytical as the derivative of the argument of tangent function is proportional 



to the overall factor i.e., = —8^yh{l — h). Finally 

1 

which is finite only for 



F(l) = In cos (114) 



1 , , 

cosK7r>0 =^ '^^2 v-'--'-^) 

Hence, for k > (1/2) the flux is not well defined on the whole segment h G [0, 1] which, in 
consequence, excludes such solutions. 

For K = 1/2, we get that F{1) = oo which leads to flux quantization and is not compatible 
with the boundary conditions for a soliton, as was discussed in Section IV.D. Solitons may, 
therefore, only exist for k < (1/2). Let us notice that for solutions corresponding to the 
choice of the model parameters approaching the critical value k = 1/2 we asymptotically 
(with arbitrary accuracy) tend to the quantized magnetic flux. The approach is rather weak 
(quadratical) as 

Ooo = — 1 + cos Kvr (116) 

which should be contrasted with typical exponential (or even more rapid) approach for the 
large gauge coupling limit in the gauged (Maxwell) BPS baby Skyrme model. Interestingly 
enough, the asymptotic quantization of the flux occurs for finite values of the parameters 
of the model which again differs from the BPS baby Skyrme model with the Maxwell field, 
where such a limit is realized only for infinite value of the gauge coupling constant. 
We remark that for generic potentials, the occurrence of zeros of the superpotential deriva- 
tive, Wh = 0, at some point Hq G (0, 1) implies that W^h becomes singular at this point, and 
the superpotential equation cannot be integrated beyond this point. In these cases, BPS 
solitons obviously do not exist. Only for some specific potentials (like the one studied in this 
section), the superpotential exists in the whole interval even if the interval contains points 
where Wh = 0, and the possible existence of BPS solitons requires some further study. 



H. Numerical results 

In a next step, we want to numerically integrate the field equations for the axially sym- 
metric ansatz, where we integrate directly the system of three equations fl69|) . fITOl) and flTTj) 
for h, a and b. Concretely, we integrated the system for the following three potentials. 
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1. The old potential V = 2h 

Solutions are compactons, and we perform a shooting from the compacton boundary. For 
the correct boundary conditions at the boundary xq, the expansions read 



,2 



a(a;) ~ ao + 0(xo - x)^, (118) 

(the coefficient /i2 here is a factor of 2 bigger than in Eq. (185!) . because here = 2/i, whereas 
we assumed V ^ h in Section III.D). The free constants are and Xq, and they must be 
determined such that the two conditions h{0) = 1 and a(0) = are satisfied. The result of 
the numerical integration for a specific choice of the coupling constants is shown in figure [H 
The values of the free constants and the energy for this solution are 

ao = -0.376724, xq = 1.5089, Eo = 15.1067. (120) 

Further, the energy coincides exactly with the BPS energy, and is exactly equal to the 
value given by Eq. (l92l) . which demonstrates that this solution is, in fact, a BPS solution. 



2. The new potential V = h{l — h) 

Solutions are again compactons, so the above discussion just repeats is now a factor 
of 2 smaller, i.e., exactly equal to (1851) . because now V ^ h). The result of a numerical 
integration is shown in figure O The free parameters and the energy are 

ao = -0.555984, xq = 4.439, Eq = 5.62986, (121) 

where ao and Eq again coincide with their BPS values. 



3. The potential V = h'^ 

Here, the Skyrme field approaches its vacuum value exponentially, therefore we perform 
the shooting from the center. The expansion at the center gives 



h{x) ~ 1 + hix + . . . 

2760/ii 2 , 
a{x) ~ X + . . . 
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h, h' a, a' 




(a) Function h with its derivative. (b) Function a and its derivative. 



5 r 



b, f 
1.5 r 




(c) Function b and its derivative. (d) Energy density. 



FIG. 1: Static solution for the old baby potential, with the coupling constant values 

/i = 7 = A = p = 'n = l. 



h{x) ~ 6o -X + . . . 

P 

where the free constants are hi and h^. For the numerical solution shown in figure [3] these 
constants and the energy take the values 

hi = -0.7071, bo = 0.663114, Eq = 8.333 

and the energy is, again, equal to its BPS value. So, the explicit numerical calculations 
completely confirm our analytical results. 



V. CONCLUSIONS 



In this paper we investigated in detail some Lifshitz field theories, characterized by an 
anisotropic scaling between space and time, with dynamical critical exponent z = d, where 
d is the dimension of space. We explicitly considered the case d = 2 for reasons of sim- 
plicity. More concretely, we studied Lifshitz field theories which have a large (infinite- 
dimensional) group of base space symmetries, namely the Special Diffeomorphisms on d- 
dimensional space. One reason to consider these theories is that the SDiff group forms a 
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relevant (infinite-dimensional) subgroup of the symmetry group of Horava-Lifshitz gravity 
0], so these models may hopefully shed some light on non-perturbative properties of the 
latter theory. Another reason is that the presence of such a large group of symmetries 
typically allows to find exact solutions of the field equations, providing thereby analytical 
nonperturbative informations about these theories. For the class of theories considered, the 
relevant term in the lagrangian density (corresponding to a renormalizable interaction w.r.t. 
Lifshitz scaling) could be interpreted as the square of a topological charge density, so that 
topological soliton solutions can be expected. Concretely, we investigated theories where the 
fields had the topology of skyrmions, but, as said, a generalization of our results to the cases 
of vortices or monopoles should pose no problem. We found that these theories have a BPS 
bound and topological solitons, i.e., skyrmions, saturating the bound. The static sector of 
;he theories is, in fact, equivalent to the static sector of the so-called BPS Skyrme models 



ul- 



161 ]. where BPS skyrmions are known to exist. Time-dependent (e.g. Q-ball) solutions. 



on the other hand, behave qualitatively different in the Lorentz-invariant BPS Skyrme mod- 
els and in the Lifshitz field theories considered here. Finally, we studied the case where the 
Lifshitz field theory is coupled to an abelian gauge field, where both scaling symmetry and 
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h, h' a, a' 




(a) Function h with its derivative. (b) Function a and its derivative. 



1.5 r 




(c) Function b and its derivative. (d) Energy density. 

FIG. 3: Static solution for the potential V = h"^, with the coupling constant values 

fi = 'y = X = p = n = l. 



renormalizability required a Chern-Simons instead of a Maxwell gauge theory. We found 
that the resulting theory, again, has both a BPS bound and solutions saturating this bound 
where, as in the case of the gauged BPS baby Skyrme model [itI, the derivation of the BPS 
bound required the introduction of a superpotential, as in supergravity coupled to a scalar 
field 
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The strongest motivation for the recent interest in Lifshitz field theories is probably 
related to the possiblitiy to find a renormalizable and unitary quantum theory of gravity, as 
proposed in Horava-Lifshitz gravity 0]. But if an anisotropic scaling is assumed for gravity 
in the UV region, then most likely the universal character of gravity makes it necessary 
to extend this anisotropic scaling also to the remaining fields and interactions. This then 
requires a detailed study of the perturbative and nonperturbative aspects of the resulting 
Lifshitz field theories, both for a check of the consistency of the whole construction and in 
order to understand whether and to which degree Lifshitz field theories may give rise to new 
physical phenomena not present in the more standard theories with invariance under Lorentz 
transformations. The results of the present paper are part of this ongoing investigation. 

The are several obvious directions in which the present work can be continued. First of 
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all, as the Lifsliitz BPS baby model has a standard time evolution, one may try to analyze 
the dynamical properties of the model. Since the model defines a well-defined Cauchy 
problem, one may, e.g., integrate the time evolution using some numerical methods. For 
example, one may be interested in scattering processes of the (exactly known) solitons or Q- 
balls. Such an analysis becomes even more relevant taking into account the infinitely large 
symmetry group of the model. Indeed, we have infinitely many conservation laws which 
may strongly influence the solitonic collisions. Hence, one may hope that a link between 
the standard integrability (e.g., manifested by a lack of radiation during scattering) and the 
generalized integrability in more than 1+1 dimensions can be found. The previously known 
examples of models integrable in the generalized sense do not have (globally) well-possessed 
Cauchy problems (Aratyn-Ferreira-Zimerman model [si], BPS baby and BPS Skyrme) or 
are defined as constrained systems (integrable submodels) by some additional derivative 
dependent constraints. Therefore, it is very difficult to investigate dynamical properties of 
these theories, either analytically or numerically. As a consequence, it has not been verified 
yet whether the scattering of the corresponding solitons occurs in the same way as in the 
usual integrable systems (no radiations, the same particles in the initial and final states 
etc.). 

Another straightforward generalization is the construction and analysis of the analogous 
Lifshitz Skyrme model in (3+1) dimensions. 
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